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We present an exploratory study of a gauge-invariant non-perturbative renormalization technique. The renormal- 
ization conditions are imposed on correlation functions of composite operators in coordinate space on the lattice. 
Numerical results for bilinears obtained with overlap and 0(a)-improved Wilson fermions are presented. The 
measurement of the quark condensate is also discussed. 



Renormalization constants (RC) of lattice op- 
erators are crucial inputs to obtain physical quan- 
tities from the matrix elements of these operators. 
We have presented in this conference a numerical 
study of a non-perturbative (NP) technique pro- 
posed in pp. This method has the advantages of 
being gauge- invariant and free of contact terms, 
in addition with the very easy implementation on 
a simulation. The first two properties can be of 
vital interest to renormalize the operators which 
enter the K ^ irir decays. 

A very exploratory analysis was already pre- 
sented in the Lattice 2002 conference |2| for bilin- 
ears operators. The results presented there were 
not conclusive and claimed for a more detailed 
study, which is the purpose of this contribution. 
We have considered only bilinears operators to 
test the method, since they are the simplest op- 
erators for which there exist results from others 
NP techniques to compare with. 

The method consists on imposing renormal- 
ization conditions on the bare green functions, 
Gr{x) = (Or (a;) Op (0)), where Or{x) is a quark 



bilinear operator, Or{x) — ipi{x)T'ip2{x)- One 
can choose 

ZUf^ = l/x) Grix) = Gr(x)(/-^ ) , (1) 

in a window 1/a ^ l/x ^ ^qcd] we will call 
this definition as the X-space scheme. One can 
also obtain Zr(M — ^/x) directly in another con- 
tinuum scheme, such MS, by imposing, 



Zj^is (^ = l/a;)Gr(x) = Gr(/i,x)(^s) 



(2) 



* Presented by J. Reyes at "Lattice 2003", Tsukuba. 
tThis work has been supported in part by the EU IHP 
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We have analyzed results for 0(a) improved Wil- 
son fermions at /? = 6.45 [a'^ = 3.87(19) GeV) 
on a 32'^ x 70 lattice and overlap fermions at 
(3 = 6.00 (a-i = 2.00(10) GeV) on a 16^ x 32 
volume. All the results presented here have been 
extrapolated to the chiral limit (further informa- 
tion can be found in [S]). 

As an example of the general strategy adopted, 
we explain the analysis for the vector-vector cor- 
relator with Wilson fermions. In Fig. ^) we show 
the dependence of the correlator (Oy(x)Oy(O)) 
on . The naive dimensional behaviour [x"^ ) 
is well satisfied but the data points are quite dis- 
persed (see also Fig. ^d)). To clarify the ori- 
gin of these effects we show in Fig. ^) the ra- 
tio between the free correlator on the lattice and 
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Figure 1. 

in the continuum, both in infinite volume (note 
that in this case the expected value is 1 up to 
0{a) effects). Comparing both figures we observe 
that the tree level calculation manifests the same 
pattern that the interacting case. These consid- 
erations together with similar analysis (see 
brought us to conclude that we are observing dis- 
cretization effects, which we propose to reduce 
considering the following ratio: 



zl{^i = l|x) = 



Gr{x) 



Gt{x) 



{F. lat.) 



Gr(x)(^-) Gr(a;)(^-^°"-) 



(3) 



where, I., F. mean interacting and free respec- 
tively. In this way we obtain the result plotted 
in Fig. which can be compared directly with 
Fig. ^p). It is clear that in this way a great part 
of 0{a) effects are canceled out in the ratio of eq. 
(PJ- The same procedure has been carried out for 
the other bilinears and for the overlap data. 

As commented above one can also obtain Zq 
directly in another continuum scheme (MS), by 
imposing (0). Now one can evolve the RCs with 
the MS Renormalization Group (RG) factors, up 
to a fixed scale, say 2 GeV. 



Z^S(2GeV) 



Z{x) 



W{^i,2GeY)R{pi,l/x) 



(4) 



where -R(^, l/x) is the matching between the X- 
space scheme and MS, which it has been calcu- 
lated at NLO H- W[^J^, 2 GeV) is the RGB factor. 
In principle /i is arbitrary as far as /i in or- 

der to avoid large logarithmic contributions. This 
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Figure 2. 



scale dependence should cancel in a RGI quan- 
tity up to truncation of the perturbative series. 
In practice we have chosen fj, — 1/x and ^ — 2/x 
and included the difference as a systematic error 
due to NNLO effects. 

Our results are summarized in Fig. [3 for the 
Wilson |(5Ji) and overlap ^p) simulation, respec- 
tively. For the Wilson case we show, as an ex- 
ample, the RCs for the scalar density and vector 
current. We also show the ratio ZpjZg. This 
ratio (as well as Zy IZa) allows to extract infor- 
mation on the renormalization window, since it 
should be a plateau up to 0(a) and NP effects. 
From Fig. [2^) we see that in the range x^ — 9-^21 
we observe a reasonable good plateau: For a;^ < 9 
large 0{a) effects are present, and above x^ ~ 21 
the presence of NP effects begins to be visible (see 
Zp/Zs)- Therefore, we choose to obtain our final 
results for the RCs in the above range. We also 
plot the scalar and vector RCs, the first in MS 
at 2 GeV and in the X-space scheme. For Zy we 
see a good plateau as one would expect. ^5(1/2;) 
in the X-space exhibits an anomalous dimension 
behaviour as it should be. On the contrary in the 
value of ^^^(2 GeV) one expects that the anoma- 
lous dimension cancels (up to NNLO correction) 
and indeed a plateau is observed. 

For the overlap simulation, we show the cur- 
rents RCs, together with the ratio Zv/Za- Since 
chirality holds exactly, one expect Zy = Z a up to 
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Table 1 

Wilson results and comparison with other methods: RI-MOM [5| and Schrodinger Functional 

Zy Z\r Za Zp/Zs Zg Zp Zt 

This work 0.801(2)(18)(6) 0.833(2)(27)(6) 0.888(2)(8) 0.702(4)(27)(23) 0.624(3)(19)(21) 0.895(2)(21) 
RI-MOM 0.803(3) 0.833(3) 0.897(4) 0.679(8) 0.609(8) 0.898(6) 

SF 0.808(1) 0.825(8) 0.912(9) — 0.61(1) 



NP effects. From Fig. I^t) one see that indeed this 
is the case for < 10, and we choose the range 
6 10 to obtain the values of Zy, Za- Contrary 
to the Wilson results, no renormalization window 
has been found for the other bilinears in the over- 
lap case. This is due to the larger lattice spacing 
(twice what we have for Wilson data) which con- 
sequently makes much narrower the renormaliza- 
tion window. In Table 1 we present our results for 
the Wilson case and those with other NP meth- 
ods. AU the RCs are in MS at 2 GeV. The first 
error in our results is statistic, the second one is 
an estimate of discretization errors and the last 
one comes from NNLO effects due to the choice 
of the scale (see above). For the overlap case we 
quote the results obtained for Zy and Za'- 

Zv ^ 1.471(2)(10)(20) Za = 1.494(2)(9)(20), (5) 

the tiny difference between the results for Zy and 
Za is due to small NP effects. In any case they are 
compatible within the errors. Our results for Za 
is comparing with Za = 1.55(4) obtained with 
the WI method in |^. From Table 1 and the 
above results in the overlap case we conclude that 
although the errors are, in most cases, larger than 
those obtained from other NP techniques, good 
agreement is found for all the cases which we can 
compare. We stress that one could reduce the 
errors earring out a recursive volume matching 
technique. 

We present a preliminary result for the chiral 
quark condensate, obtained in a novel way. The 
idea is very simple and consists on studying the 
chiral limit of the trace of the quark propagator 
in X-space together with its OPE: 



(6) 



We stress that in once the l.h.s. is renor- 
malized, the r.h.s. gets renormalized too, and no 
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Figure 3. ^ 

other divergence (apart from the natural one at 
a; = 0) is present, thus only the renormalization 
constant for the quark field is needed. We have 
computed the needed Wilson coefficients factors 
in eq. ©. From eq. ((BJ and in the chiral limit 
one can obtain an estimate of the chiral conden- 
sate. This is shown in Fig. 13 From that figure we 
see large discretization effects for < 15. One 
can see also that there is a contribution from the 
next order term which goes as - From a linear 
fit we obtain our (preliminary) result: 



(VV') "(aCeV) = -(276(5)(13)MeV)' 



(7) 



in good agreement with other lattice determina- 
tions. The first error is statistic and the second 
is due to the uncertainty in the value of a^^. 
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